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Field-Induced Quasiparticle Excitation in Ca(Al sSio.s^: Evidence for 
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The temperature (T) and magnetic field (H) dependence of the magnetic penetration depth, 
X(T,H), in Ca(Alo.5Sio.5)2 exhibits significant deviation from that expected for conventional 
BCS superconductors. In particular, it is inferred from a field dependence of \(H) (oc H) at 
2.0 K that the quasiparticle excitation is strongly enhanced by the Doppler shift. This suggests 
that the superconducting order parameter in Ca(Alo.5Sio.s)2 is characterized by a small energy 
scale As/ks < 2 K originating either from anisotropy or multi-gap structure. 
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Since the discovery of binary inter/metallic compound 
superconductor MgB 2 (T c =39 K), 1 the origin of pair cor- 
relation leading to such high critical temperature has 
been drawing much interest. Up to now, a majority of 
experiments as well as theories suggests that the pair cor- 
relation is mediated by a strong electron-phonon interac- 
tion, where the relatively high T c is explained by the light 
mass of the two-dimensional honeycomb layer formed by 
B atoms. The two-dimensional feature is also preferable 
for the strong electron-phonon interaction, and thereby 
the AlB2-type crystal structure common to MgB2 is cur- 
rently attracting attention as a possible basis for devel- 
oping new superconductors. 

In binary silicides, ThSi2, 2 USi2, 3 and several rare- 
earth metal disilicides have the AlB2-type structure. j3- 
ThSi2 is known to be a superconductor with a crit- 
ical temperature T c =2.41 K. 4 Above 16 GPa, CaSi 2 
takes a AlB2-like structure, and it becomes supercon- 
ducting with T c =14 K. 5 Recently, a new ternary silicide, 
Sr(Ga,Si)2, which also has the AlB2-type structure, was 
reported to be a superconductor with T c =3.4 K, 6 stim- 
ulating active investigation of analogous compounds. In 
this class of materials, ternary silicide Ca(Aln. 5810.5)2 has 
the highest critical temperature T c =7.7 K. 7 It is reported 
on this compound that the behavior of electron-heat ca- 
pacity deviates from that of the BCS-type, and that the 
effect of hydrostatic pressure on T c is positive. 8 Unfortu- 
nately, despite various experiments so far, there is very 
little known on the structure of superconducting order 
parameter in Ca(Aln. 5810.5)2- 

In this paper, we report the quasiparticle excitation in 
the flux line lattice (FLL) state of Ca(Aln. 5810.5)2 studied 
by muon spin rotation (/iSR) . The magnetic penetration 
depth A, which reflects the quasiparticle excitation, can 
be determined by /iSR due to the inhomogeneity of the 
magnetic field distribution in the FLL state. We show the 
temperature (T) and magnetic field (H) dependence of 
A in polycrystalline sample of Ca(Al . 5810.5)2- Our result 
indicates that A is proportional to T 3 at low temperature, 
and that it exhibits a steep increase with almost linear 



dependence on the applied magnetic field. These findings 
strongly suggest that Ca(Al . 5810.5)2 has an anisotropic 
order parameter or a multi-gap structure with a band 
having a small gap energy, which is qualitatively similar 
to the case in MgB2. 

The samples were prepared by arc melting method 
with a stoichiometric mixture of Ca (99.9%), Al 
(99.99%), and Si (99.999%) (in 1:1:1 composition) un- 
der an Ar atmosphere. The structure of polycrystalline 
sample was examined using a powder x-ray diffraction 
technique. Diffraction signal from Ca(Aln. 5810.5)2 phase 
was observed as a main contribution besides those from 
a small amount of impurity phases. Magnetic suscepti- 
bility and electrical resistivity measurements were per- 
formed with the SQUID magnetometer (MPMSR2) and 
the PPMS system (Quantum Design Co. Ltd.). The 
residual resistivity ratio was estimated to be 2.85. 

Figure 1 shows the temperature dependence of mag- 
netic susceptibility at 10 Oe. The solid triangles and cir- 
cles show data obtained in zero-field cooling (ZFC) and in 
field cooling (FC), respectively. In FC, a Meissner effect 
can be seen below 7.7 K. The inset of Figure 1 shows the 
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Fig. 1. Temperature dependence of magnetic susceptibility at 10 
Oe. The inset shows the upper critical field versus temperature. 
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temperature dependence of the upper critical field (H c2 ) 
which was estimated from the temperature dependence 
of electrical resistivity at respective magnetic hclds. The 
solid line is the result of fitting analysis by the following 
equation from a local-paring theory, 



H c2 (T) = H c2 (0)(l 



_ t 3/2n3/2 



(1) 



where r = T/T c , which yields H c2 (0) = 2.60(2) T. 

The /zSR experiment was performed on the M15 beam 
line at the Tri-University Meson Facility (TRIUMF, 
Canada) which provides a muon beam with a momentum 
of 29MeV/c. The specimen having a dimension of about 
7x7 mm 2 was mounted on sample holder and placed in a 
cryostat. The sample was field-cooled at every magnetic 
field points to minimize the effect of flux pinning. The 
temperature and magnetic field dependence of transverse 
field (TF-) ^SR spectrum was obtained at H = 0.05 T 
and T = 2.0 K (=0.26T C ), respectively. Since the muon 
stops randomly on the length scale of the FLL, the muon 
spin precession signal P(t) provides a random sampling 
of the internal field distribution B(f) 



P(t)=P X {t)+lPy{t) 



i{B) exp{i lfl Bt+c/>)dB, (2) 



S 
< 




-120 -80 -40 40 80 120 

B-H [G] 

Fig. 2. Fast Fourier transform of /^SR spectra in Ca(Alo.5Sio.s)2 
at 2.0 K under several applied magnetic fields. 



n(B) = (S(B(f) - B)) r , 



(3) 



where 7^ is the muon gyromagnetic ratio (= 2n 135.53 
MHz/T), (f> is the initial phase of rotation, and n(B) is 
the spectral density for the muon spin precession sig- 
nal determined by the local magnetic field distribution. 
These equations indicate that the real amplitude of the 
Fourier transformed muon spin precession signal corre- 
sponds to the local field distribution n(B). The London 
penetration depth in the FLL state is related to the sec- 
ond moment of the field distribution ((AB) 2 }, which is 
reflected in the ,uSR line shape. 9 For polycrystalline sam- 
ples, the Gaussian distribution of local field is a good 
approximation, 



P(t) ~ exp(-a 2 r/2)exp(i7^iJt), 



(4) 
(5) 



For the ideal triangular FLL with isotropic effective car- 
rier mass to*, A is given by the following relation, 9-11 

a^s' 1 ] = 4.83 x 10 4 (1 -h)[l + 3.9(1 - fc) 2 ] 1/2 A- 2 [nm]. 

(6) 

where h = H/H c2 (T), and A is given by the following 
relation of superconductive carrier density n s (T, H) 



A 



TO C 



(7) 



4>im a (T, H)e 2 

It should be noted that the reduction of a with increasing 
field described by Eq. (6) is due to the stronger overlap 
of B(f) around vortices at higher fields while A is a con- 
stant. The reduction of n s (T, H) due to the quasiparticle 
excitation leads to the enhancement of A, which causes 
further reduction of a from that expected from Eq. (6). 

Figure 2 shows the fast Fourier transforms (FFT) spec- 
trum of muon spin precession signal in Ca(Alo.sSio.5)2 
for several magnetic fields at 2.0 K, where the real ampli- 
tude of FFT corresponds to the internal field distribution 



n(B) in the FLL state convoluted with additional damp- 
ing due to small random field from nuclear moments. As 
most explicit in FFT spectrum at 0.05 T, the line shape 
is characterized by the shift of peak to a lower field and 
also associated broadening of linewidth due to the for- 
mation of FLL. 

Considering possible contribution from muons stop- 
ping in the normal part of the specimen, we adopted 
two components with the following empirical form in an- 
alyzing the /iSR time spectra, 



AP(t) 



2 

A i ex P ( - 



2 t 2 



exp( 7AtJ ff l t + fa) (8) 



where the index i refers to the components of supercon- 
ducting (i = 1) and normal (i = 2) domains, A is the 
total positron decay asymmetry with Ai being the par- 
tial asymmetry, <7j is the muon spin relaxation rate, Hi 
is the central frequency, <pi is the initial phase for respec- 
tive components. The superconducting volume fraction 
(= Ax/ (At + A 2 ) was obtained to be ~ 0.95. 

The temperature dependence of o\ (oc n s ) is shown 
in Figure 3. In the region below T c where the FLL is 
formed, o\ increases with decreasing temperature. Ac- 
cording to the empirical two-fluid model (which is a good 
approximation of the BCS theory), the following relation 
is expected to hold, 

CTl (T) = CTl (0)[l - r 4 ] cx n s (T) cx A" 2 (T). (9) 

Solid line is the result of fitting analysis by a similar 
formula with an arbitrary power, 

<7 1 (T)=<7 1 (0)[1-^], (10) 

wheh yields (i = 2.73 and T c = 7.53(1) K when both 
P and T c are assumed to be free parameters. The result 
also means that the deviation of A, 



<5A = A(T) — A(0), 



(11) 
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Fig. 3. Temperature dependence of the Gaussian relaxation rate 
<ti at 0.05 T. The solid curve is a fitting result with cr\(T) = 
ct(0)[1 — (T/Tc)' 3 ], while the dashed curve corresponds to the 
case /3 = 4. Inset: A plotted against T 3 . 



exhibits a tendency predicted for the case of line nodes 
(d-wave pairing) with some disorder (i.e., dirty limit, 
where S\ oc T 2 ). 12 Since the gap energy A(T) ~ A (l - 
t 4 ) is predicted to be least dependent on temperature 
for r < 0.4 in the BCS model, the observed deviation 
of f3 from 4 suggests the presence of excess quasiparticlc 
excitation. This leads to the possibility that the order 
parameter has an anisotropic strutcure, where the ex- 
cess quasiparticles are induced at the vicinity of nodes 
(\Af,\/kB < 2 K). Another possibility may be that, as- 
suming a multi-gap structure with As being one of the 
smallest gap energy, the excess quasiparticles are due to 
the thermal activation over As (As/fee < 2 K) ; unfor- 
tunately, the absence of data below 2.0 K did not allow 
us to perform the reliable fitting analysis using a two-gap 
model. 15 

The presence of an energy scale smaller than that 
for the single BCS-gap energy Ao is further supported 
by the magnetic field dependence of A. As shown in 
Fig. 4a, a(h) decreases with increasing field much steeply 
than expected by Eq. (6). Accordingly, A exhibits a 
strong field dependence, where A increases almost lin- 
early with h (see Fig. 4b). This is similar to the case of 
NbSe2 16 and YNi 2 B 2 C 17 which exhibit characters spe- 
cific to anisotropic order parameters in spite of the sug- 
gested s-wave symmetry, or to that of MgB 2 18 having a 
two-gap structure. Solid curves in Fig. 4 is the result of 
fitting by the following linear relation, 



X(h) = A(0) [1+77/1], 



(12) 



where 77 is a dimension-less parameter which represents 
the strength of pair breaking effect. Fitting yields 77 = 
0.80(8) (with H c2 = 2.12(2) T determined by Eq. (1)) 
which is comparable with that in NbSe2, YNi 2 B 2 C and 
MgB 2 ( i.e., ry = 1.61, 0.97, and 1.27, respectively). 

In the FLL state, the quasiparticle energy spectrum 
around the vortex cores is shifted by an amount e = p ■ v s 
due to the semiclassical Doppler shift, where p is the 
quasiparticlc momentum (~ 771*7)/?), and v s is the ve- 
locity of supercurrent. Since the density of state (DOS), 
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Fig. 4. Magnetic field dependence of the Gaussian relaxation rate 
cri (a), and penetration depth A (b) at 2.0 K. Solid curves are 
fitting result with \(h) = A(0)[1 + r/h], and the dashed curves 
correspond to r\ = 0. 



N(E), for anisotropic order parameters takes a non-zero 
value for an energy E > , quasiparticles can be ex- 
cited by the Doppler shift outside of the vortex cores 
with a population proportional to N(e), leading to the 
enhancement of A. 19 The magnitude of 77 represents the 
degree for the increase of DOS for quasiparticles, which 
must be roughly proportional to the phase volume of the 
Fermi surface where the Doppler shift exceeds the gap 
energy (s > \ Ag | ) . It also follows that the effect depends 
on the direction of v s (and hence that of the external 
field H relative to the order parameter) in a single crys- 
talline specimen. According to Volovik, the quasiparticle 
density of state for anisotropic order parameter is 

A dcloc (0) ~ N F K& L Vh = N F g(h), (13) 



K (x 



|Ad<e 



lArJdk, 



(14) 



where Np is the DOS for the normal state and K is a 
constant on the order of unity. 19 It is important to note 
that K is proportional to the phase volume of the low 
excitation energy in A? , thereby carrying information on 
the degree of anisotropy for A?; the factor ft, 1 / 2 comes 
from the inter-vortex distance (oc ft -1 / 2 ) multiplied by 
the number of vortices (oc h). The superfluid density at 
a given field is then 



n s (h) ~n s (0)[l- 5 (/i)], 



(15) 



which is directly reflected in the magnetic penetration 
depth through Eq. (7). Therefore, as a mean approxima- 
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tion, we have 
\(h) = 



A(0) 



X(0)[1 + cKe GL h], (16) 



y/ 1 - 9(h) 

where c ~ 1.5 for < h < 0.5. 20 Thus, the comparison 
between Eqs. (12) and (16) yields 



(17) 



indicating that the slope r\ reflects the phase volume 
of the Fermi surface where |Ar| < e. The situation in 
the case of multi-gap state is similar when As < fc^T. 
For example in the two-gap model, 15 using the relative 
weight, 1 — x, of the quasiparticle DOS for the smaller 
gap As , we have 



Kotl 



x. 



(18) 



Since the Doppler shift is far smaller than the gap en- 
ergy in usual situation for the isotropic gap, no such en- 
hancement is expected for the conventional s-wave pair- 
ing (i.e., 77 <Cl for |A |, |A S | » k B T). For example 
in Y(Ni .8Pto. 2)2620, which behaves as a conventional 
BCS superconductor with an isotropic gap, it is reported 
that r] ~ 0. 21 This is in good contrast with the case 
of YNi 2 B 2 C in which an anisotropic gap (two-gap state 
with point nodes) is strongly suggested by other experi- 
ments. 22 ' 23 On the other hand, a stronger magnetic field 
dependence is predicted for the case of <i-wave pairing 
due to the relatively large value of K. Typical examples 
of the d-wave pairing are found in high-T c cuprates, in 
which 77 is reported to be 5-6.6 for YBCO. 24 The com- 
parison of those earlier results with our result suggests 
that Ca(Al .5Sio.5) 2 has an anisotropic order parameter 
comparable to that of YNi 2 B 2 C or a multi-gap order pa- 
rameter in which one of those has a small gap energy 
(I A s I As < 2 K) as found in MgB 2 . 

Unfortunately, our results were obtained using poly- 
crystalline samples, which makes it difficult to deduce 
A by analyzing data using the well-defined microscopic 
model. Early experiments on high quality single crys- 
tal sample of Ca(Al .5Sio.5) 2 have reported that the up- 
per critical field has an anomalous angular dependence 
which deviates from the Ginsburg-Landau anisotropic 
mass model. 25 Moreover, very recent report suggests that 
the crystal structure has clear fivefold and sixfold super- 
lattice. Therefore, we are preparing high quality single 
crystals of Ca(Al . 5810.5)2 for further /^SR study of this 
compound in more detail. 

In conclusion, We have performed /iSR experiments on 
polycrystalline Ca(Alo.5Sio.s) 2 to elucidate the structure 
of superconducting order parameter. We found that n s (t) 
exhibits a considerable deviation from that expected for 
the BCS prediction, and that X(h) exhibits significant in- 
crease with applied magnetic field. These results strongly 



suggest the presence of unconventional order parameter 
in Ca(Al . 5 Si .5)2. 
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